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Abstract 

Motivated by the Goldbach conjecture in Number Theory and the abelian bosonization 
mechanism on a cylindrical two-dimensional spacetime we study the reconstruction of a real 
scalar field as a product of two real fermion (so-called prime) fields whose Fourier expansion 
exclusively contains prime modes. We undertake the canonical quantization of such prime fields 
and construct the corresponding Fock space by introducing creation operators bj, -labeled by 
prime numbers p— acting on the vacuum. The analysis of our model, based on the standard rules 
of quantum field theory and the assumption of the Riemann hypothesis, allow us to prove that 
the theory is not renormalizablc. We also comment on the potential consequences of this result 
concerning the validity or breakdown of the Goldbach conjecture for large integer numbers. 
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Omnia disce, videbis postea nihil esse superfluum 
(Learn everything, later you will see that nothing is superfluous.) 
Hugh de Saint Victor (c. 1078-1141), French or German philosopher 



1 Introduction 

In this paper, we focus on the border between physics and number theory. It can be regarded as 
a first step to investigate, in particular, whether quantum field theory can shed any light on some 
long-standing, yet unproved, mathematical problems in number theory. Indeed it seems reasonable 
that the "translation" of certain difficult mathematical questions into the physical language might 
provide essential clues to better understand (or maybe help to solve) them. For instance, it is 
well-known that the Hamiltonian describing quantum-mechanically a one-dimensional harmonic 
oscillator has an equally spaced spectrum (in the appropriate units and after a shift corresponding 
to the zero point energy) given by the set of positive integers. This remarkable property naturally 
leads to the following question: Is it possible to use spectral properties of this or other operators 
to study and analyze problems of number theory, e.g. the distribution of primes or the Goldbach 
conjecture? 

At the beginning of the 20th century, Polya and Hilbert independently suggested the possibility 
of proving the Riemann hypothesis on the distribution of non-trivial zeroes of the Riemann function 
£(s) by finding a Hermitian operator with eigenvalues that could be put in correspondence with 
those of £(s) OP- Moreover, a remarkable connection between number theory and nuclear physics 
emerged (see [2] for a review) as a consequence of a fortuitous meeting at Princeton in the early's 
1970s between the mathematician Hugh Montgomery and the physicist Freeman Dyson. It was 
then noticed that Random Matrix Theory (introduced in the Physical Sciences by Eugene Wigner 
in the fifties OH]) could model apparently unrelated problems of both fields. The level spacing of 
heavy nuclei and the non-trivial zeroes of the Riemann zeta function showed a striking statistical 
similarity. 

This example is not the only one to suggest a link between physics and prime numbers. In 1990 
Bernard Julia put forward the idea of a fictitious, non-interacting, boson gas [5] whose constituents 
were called primons. The partition function of the primon gas was found to be the Riemann-zeta 
function C( s )- The idea of the Riemann gas can also be extended to fermions leading to a relatively 
simple partition function of the form £(s)/£(2s). Furthermore, it has been argued that more 
intriguing relations exist between the Riemann gas, quark confinement and even string theory [6]. 
This is a clear example of how studies initially started as a rather unorthodox approach to the 
Riemann hypothesis can lead, in turn, to new and eventually fertile perspectives in physics and 
vice versa. A thorough summary of several suggestive relationships between different branches of 
physics and the Riemann hypothesis can be found in (7j. 

In this paper we develop a (quite simple) physical model involving a free fermionic field defined 
on a cylindrical space-time manifold, with prime excitation modes only. No relation with p-adic 
theories in physics [8j[9] has been envisaged in this work. Rather, one of our goals is to introduce 
prime numbers within the standard language of Quantum Field Theory (QFT) and investigate its 
primary consequences. In particular, we will cast Goldbach's conjecture into the standard language 
of Quantum Field Theory (QFT). Interestingly, as will be shown, the consistency analysis for the 
underlying QFT ultimately involves the Riemann hypothesis in the computation of the renormalized 
vacuum energy. 
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1.1 Goldbach's conjecture 



In June 1742 the mathematician Christian Goldbach wrote a letter to his friend, Leonhard Euler, 
where he made his famous conjecture about prime numbers which can be split into the following 
two statements: 

Every even integer greater than 2 can be expressed as the sum of two primes 

Every odd integer greater than 5 can be expressed as the sum of three primes 

The first statement is known as the binary or strong Goldbach's conjecture (GC). The second is 
known as the ternary or weak GC. Note that if the strong GC is true then the weak GC is also 
true. If not otherwise stated, we will denote the strong GC as the GC for short@ 

The weak GC was proved in 1923 by Hardy and Littlewood, for sufficiently large odd numbers 
as an "asymptotic" result, assuming the truth of the Riemann hypothesis. In 1937 Vinogradov 
again proved the weak conjecture for a sufficiently large (but indeterminate) odd number using 
analytic methods. In 1966 Chen Jing-Run proved that every sufficiently large even number can 
be expressed as the sum of a prime and a number with no more than two prime factors. At the 
present moment, the GC has been checked up to integers ~ 10 18 [10] by using computers. 

1.2 Polignac's conjecture 

In 1849 Alphonse de Polignac (1849) made the general conjecture that there are infinitely many 
cases of two consecutive prime numbers such that their difference is any even number n even . For 
n>even = 2 it reduces to the twin prime conjecture. For n even = 4 there should be infinitely many 
cousin primes... Based on heuristic arguments, Hardy and Littlewood found in 1922 a law to 
estimate the density of twin primes [11] . also discussed by Clement in 1949 |12j . 

The following weakened form of Polignac's conjecture (PC) can be enunciated by relaxing the 
condition on the necessity that both primes are consecutive: 

Any even number can be expressed as the difference of two primes in infinitely many ways. 

Polignac's weakened conjecture can be seen as complementary to Goldbach's conjecture. Both 
Polignac's and Goldbach's conjectures have never been proved. 

1.3 Riemann hypothesis 

At first glance the Riemann hypothesis (likewise the two preceding conjectures) does not seem very 
complicated to state, but in fact it lies deep in the roots of number theory, in connection with the 
mysteries of the distribution of the prime numbers: 

All non-trivial zeros of the zeta function have real part equal to one-half. 

The Prime Number Theorem is equivalent to the assertion that no zeros of £(s) on the complex 
plane (s = a + it) lie on the boundary of the critical strip (0 < a < 1). Understanding how the 
non-trivial zeros are located relative to the critical line (a = 1/2), remains today one of the deepest 
unsolved questions in Mathematics. Moreover, note that many developments in Number Theory 
are contingent upon the validity of the Riemann hypothesis (RH). In this paper we explore the 
connection between the RH and the GC and PC on the light of a quantum field theoretical model 
to be developed below. 

4 The number 1 was considered as prime in Goldbach's time, but presently it is usually not. This issue is, in any 
case, irrelevant for our study. 
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2 Scalar and fermion fields on a cylindrical space-time 



Our starting point is the Lagrangian density of a classical free real scalar field in a two- 

dimensional spacetime: 



C 



(1) 



where x denotes the space coordinate, leading to the well-known massless Klein-Gordon equation: 

Q2 Q2 1 



dt 2 dx 2 



$(i,x) = 



(2) 



The inverse metric tensor g^ v has been defined such that g 00 = 1, g 11 = —1, all other components 
vanishing. In addition we will set H = c = 1 throughout the paper. Since we have assumed a zero 
mass the theory inherits the two-dimensional conformal symmetry [13], given by the spacetime 
transformations x^ — > with x^ = (t±x)/\/2. 

As we ultimately want to deal with integer numbers we will let the spatial dimension "curling up 
into a circle" , becoming cyclic thereby introducing a discrete spectrum in the theory. In such 1 + 1 
dimensional space, the scalar field §>(t,x) has to be periodic, satisfying the boundary condition: 



<f>{t,x) = <S>(t,x + 2irR) 



Let us now expand this field as 



$(t,x) 



1 



oo 



¥>n(*) e ' 



inx/R 



(3) 



(4) 



where <p n (t) are complex Fourier coefficients in the above expansion, corresponding to a given radius 
R. The coefficient in front of the sum is introduced for normalization purposes. Using the equation 
of motion ([2]) one can fix immediately the (normalized) positive frequency solutions 



¥>n(*) 



-lUJ n t 



\/2u n 



(5) 



where u n = \n\/R (for simplicity we neglect the zero-mode solutions). The modes labeled by 
positive values of n represent complex waves that move from left to right, while negative values 
represent left-moving modes. Obviously, owing to the linear character of the Klein-Gordon equation 
([2]), the $>(x,t) field can be split into left and right moving (chiral) fields, <&-(t,x) and $>+(t,x), 
respectively, both obeying the periodic constraint ([3]), as: 



$(t,a?) = *-(t,x) + *+(t,x) 



(6) 



in such a way that the former (latter) only contains (positive frequency) Fourier components labeled 
by positive (n > 0) and negative (n < 0) integer numbers, respectively. At the quantum level the 
Fourier expansion is re-expressed in terms of creation and annihilation operators 



$_(t - x) 
<S>+(t + x) 



1 



V2itR 



E 

n>0 

E 

n<0 



e ~i(uj„t-nx/R) , 



a 



-t 

« i(w n t-nx/R) 



-i(u n t— nx/R) 



+ 



^ e i(u) n t-nx/R) 



y/2 



(7) 

(8) 
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where [a~,am] = S nm , [a~,a~] = = [an\cim], and also [4,oin] = Km, [a£,a m ] = 
[a^,am]. The vacuum state |0) is defined by the relations a£\0) = 0. 

We go further and generalize the above splitting in terms of even and odd integer numbers 



$(t,x) = $ {even \t,x) + ^ odd \t,x) 



with 



^ even \t,x) 
^ odd \t,x) 



1 



V2irR 
1 



E 



a n _ e -i(ui„t-nx/R) , 



a 



t 

i(u} n t-nx/R) 



n odd 



2 e -i(unt-nx/R) 



+ 



a rc i(u n t-nx/R) 



y 7 ^ 



(9) 

(10) 
(11) 



with analogous commutation relations for the creation and annihilation operators. From the point 
of view of CFT the decomposition ([9]) is somewhat bizarre, since it partially breaks conformal 
invariance. Only the Virasoro generators L n = 1/2 V , : a m a n - m :, with n being an even 

— '"'even/odd 

integer, do preserve the above decomposition. However, in our considerations conformal symmetry 
will play a rather heuristic role. 

Now, with the goal of restating the GC and the PC in a field-theoretic language, we write the 
even (positive) numbers as either the sum or the difference of two (non necessarily consecutive) 
primes p and q: 



n. 



p + q 
\p - q\ 



(12) 
(13) 



From Eq. (|10p and taking into account the elementary property of the exponential function, one 



could make the tentative hypothesis that the chiral field $ 



{even) 



(t,x) 



(even) 



(t,x) 



V2ttR 



E 

>o 



-in(t—x)/R 



+ 



J,n(t—x)/R 



<*>r> 



(14) 



can be reconstructed as a composite field according to &^ ven \t — x) oc (f)-(t — x)cf)-(t — x) where the 
chiral field 4>- it — x) is expressed by definition as a Fourier expansion over prime modes exclusively. 
This ansatz does not fit smoothly with the underlying conformal symmetry. Neither $_ and (/>_ have 
well-defined properties under conformal transformations. However, this proposal can be improved 
on the basis of the classical formulae on the abelian bosonization in two-dimensions |14|I15|. It is 
well-known that the light-cone derivatives of the scalar field <9±$ are conformal fields which can be 
reconstructed in terms of fermion bilinears, as we review next (see also |16j). 

2.1 Bosonization in two dimensions 

To describe massless fermions in two dimensions it is useful to use a real representation of 7-matrices 



1 

and write a Dirac field using two real components 



-1 

1 



7 s = 7 <y 



1 
-1 



(15) 



(16) 
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The equation of motion = implies immediately that \I/ _ = *S>-(t — x) and = ^f + (t + x). 

Moreover, the theory turns out to be invariant under conformal transformations — > y ± (x =l= ), 
where the chiral Fermi fields transform as follows 



dx- 



1/2 



*±(y ± ) = ( =± J *±(y ± (x ± )) . (17) 

This means that the conformal algebra for the Dirac field has central charge c = 1/2. It is possible 
to reconstruct a scalar field which has central charge c = 1, from two real Dirac fields with 
i = 1,2. The currents ^Jj.^^ acquire a simple form in terms of derivatives of <5 

: ;= T -^=d±1> , (18) 

V Z7T 



where one introduces normal ordering to avoid divergences. We will also assume for simplicity that 
the fermionic fields satisfy the same periodic conditions H of ([3]), i.e. VP^_(i, a; + 2itR) = ^f-^_(t,x). 
Expanding the fermion fields in modes (for simplicity we restrict to the right-moving sector) 

* j (t-x) = -=L= V V n e -i"(fi-*)/R (19) 

where {b l n ,lP m } = <5 JJ <5 ni _ m , one recovers the scalar annihilation and creation operators as bilinears 
involving all fermion operators b % T (we use the usual definitions bp = b % _ r and = aZ n ) 

a~ = sgn(n) ^ % : b\b 2 n _ r : , (20) 

r 

with the commutation relations [a~,am] = 5 n)Tn , with n > 0, or, equivalently 



(21) 



2.2 The even scalar field in terms of "odd" fermionic fields 



The trivial fact that any even positive number can be written as a sum of two odd numbers suggests 
the possibility of extending the reduction of an even boson field in terms of two odd fermion fields. 
In fact, one can easily obtain a relation (and a similar one for the opposite chirality) equivalent to 
Eq. (|20p this time for entire n such that: 



sgnyn 



Y,i--b l r b 2 n _ r : (22) 



r dd 



where bi in the above expression are annihilation and creation operators associated with fermion 
fields with Fourier expansion over odd modes exclusively 

* {odd)j = -=L= Y il e -Mt-x)/R _ (23) 

IL odd 

The difference with respect to the standard fields is that the commutation relations for the bosonic 
creation/annihitation operators are now of the form (with n and m even entires): 



5 Fermionic fields can also satisfy anti-periodic conditions yielding half-integer modes in the Fourier expansion. 
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2.3 The even scalar field in terms of "prime" fermionic fields 

The above mechanism and the GC (and also PC) suggest to continue the reduction process and 
consider fermion fields i/P whose Fourier expansion exclusively contains (odd) prime mode Jfl 

P 

The (even) bosonic field can be reconstructed from the expression 

«n = sgn(n) ^ i : b\b 2 n _ r : , (26) 

±reP' 

and it has now the following non-canonical commutation relational 

r(n) 



n 



(27) 



with n and m even entires; r(n) stands for the number of times the positive even integer n can be 
written as the sum of two primes (the so-called Goldbach's partition). 

Notice that the Goldbach conjeture has been invoked in order to reconstruct the field Q_Q even 
as a product of two fields ^ with only prime Fourier components. This becomes explicit in the 
assumed non- vanishing function r(n) in the above commutation relations. 



2.4 Even and odd states as "prime" excitations of the vacuum 

According to quantum field theory, b P and bp = b{^ p can be interpreted as annihilation and creation 
operators acting on the vacuum state denoted as \Qr) such that 

V p \0 R ) = (28) 

Once the creation and annihilation operators have been defined, one can readily obtain two-prime 
states as 

\ Pl , P2 ) = N bl pi bl p2 |0 fl > (29) 

where n even = P1+P2', and N stands for the normalization of states: N = 1 if p\ ^ p 2 , or N = l/y/2 
if pi = p 2 - 

Similarly, a three-particle state results from 

\pi,P2,Ps) = N bl pi bP_ p2 b k _ p3 \0 R ) (30) 

with N = 1 for Pl ^ p 2 ^ P3, ■ ■ ■ satisfying n odd = p\ + p 2 + P3- 

The validity of the GC implies that every even state, whose "mass" is given by M n = n even /R, 
can be interpreted as a two-particle (two-prime) state \pi,P2), where n even = p\ +p 2 - Correspond- 
ingly, every odd state, with mass M nodd = n dd/R can be interpreted as either a one-prime state (for 
n-odd a prime number itself) or a three-prime state \pi,P2,P3), where n even = p\ +p2 +PsH- Leaving 
aside the fact that we are dealing with a noninteracting field theory, states such as the two- and 

6 From our point of view it is natural to consider the number one as a prime number. We denote by P' the set of 
all (odd) prime numbers, including the number one. 

7 One can recover standard commutation relations with a redefinition of the operators a„ ■ 

8 Actually, primes greater than 5 can be expressed as the sum of three primes according to the ternary GC. 
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three-prime excitations of the vacuum bear some resemblance to hadrons, i.e. mesons and baryons, 
interpreted in non-relativistic models of particle physics as quark-antiquark and three-quark (or 
antiquark) states, respectively. 

Naturally, four-prime states (or higher) can be formed by acting four times with the creation 
operation on the vacuum: 

\pi,P2,P3,Pa) = N bl pi bi p2 b k _ m b l _ p4 \Q R ) (31) 

and so on for any even combination of creation operators. 

Given an even integer n even in our physical model, the GC implies that all those Fock states 
such that n even = pi + P2 = v'\ + P2 + P3 + Pa = "'■> are m ^ ac * degenerate. (Of course, similar 
extensions can be applied to odd combinations of creation operators.) Would the GC fail for a 
given n even , the corresponding state could not be represented as a two-prime state anymore but 
likely as a higher Fock state. In Particle Physics jargon, one would refer to such states as "exotic". 

The particle interpretation of the model presented above requires a consistent definition of the 
vacuum energy on a two-cylindrical spacetime. We examine in detail this issue in the following 
sections. Specifically, we will study the regularization and renormalizability properties of both a 
bosonic field and a fermionic field on the circle, the latter with vanishing Fourier coefficients for 
modes with non-prime labels. We will pay special attention to the two point functions and the 
expectation value of the Hamiltonian operator (Casimir energy). 



2.5 Vacuum zero-point energy 

The expectation value of the the Hamiltonian H§ of the free scalar field 3>(cc,i), can be written as 

(0r\H*\0 r ) = - J2 w n, w n = \n\/R (32) 

n 

where the sum runs over all integer numbers and the vacuum state corresponds to a given radius 
R. This quantity obviously diverges as n — > 00. 

On the other hand, in the absence of gravity the zero-point energy actually has no physical 
meaning in itself. This motivates the removal of the infinite vacuum energy of free Minkowski 
space by defining the Casimir energy Eq as in |17H19j : 

E c = (OrI-HsIOr) - lim (0 R \H*\0 R ) (33) 

R— >oo 

Nowadays it has become standard practice the use of a regularization method to determine the 
Casimir energy based on the Riemann Q function [20|,l21j defined as 

coo = J>~ s ( 34 ) 

n 

Analytic continuation of £(s) from the region of the complex plane such that Re(s) > 1, to the 
whole complex plane including the real negative axis (and therefore s = — 1) permits the regularized 
summation of the divergent series (|32p . Thus one easily gets the textbook result for the vacuum 
energy density pq: 

Po = 2nR = (35) 
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Moreover, the vacuum energy energy for the chiral fields <&± is, obviously, half of the vacuum energy 
density obtained above for <I>. We also note that the Casimir energy density when the sum over 
the modes is restricted to even (odd) integers is -1/6ttR 2 (l/127ri? 2 )|. 

2.5.1 Fermionic fields 

For the massless Dirac field ^ considered before the expectation value of the Hamiltonian Hq, is 
minus the corresponding scalar field result 

(0 R \Hy\0 R ) = ~-^2w n , w n = \n\/R (36) 

n 

Therefore, the regularized vacuum energy density turns out to be positive, with value l/24irR 2 . 

On the other hand, we are considering ip J (t,x) introduced in Eq. (|25p as the fundament fields 
in our model. Consequently, the expectation value of the hamiltonian reads now 

(0r\H^\0 r ) = -- w P , w p = p/R (37) 

pGP' 

where the sum runs over all (odd) prime numbers. Similarly to (|32p and (j36[) . this quantity also 
diverges as p — > oo and requires a regularization procedure. 

In fact, the formal similarity of summations (|32|) and (|37p suggests trying a similar regularization 
as before but this time using the the Riemann Prime zeta function, defined aJ^l 

P(s) = Y,P~ S - (38) 

Unfortunately, P(s) cannot be analytically continued to the negative real half-plane (Re(s) < 0) 
because of a natural boundary along the line Re(s) = [23]. Therefore, the above method to 
obtain the Casimir energy, though successful for integer modes as shown before, cannot be applied 
to our model with prime modes only. This failure should not be seen as exceptional as the Casimir 
effect regularization based on the use of the analytic continuation of the Riemann function is 
known to be not always applicable. This happens, for instance, in the case of a massive field on 
a two-dimensional cylinder [23]. In the next section we will approach this important question by 
a more systematic method based on the point-splitting regularization and renormalization tools of 
quantum field theory [T8 l l25 | l26] . 

3 Quantum prime field theory 
3.1 Scalar field 

It is instructive to consider the calculation of the vacuum energy of the scalar field &(t, x) by means 
of the so-called Abel summation of the divergent series ([32]) . To motivate this procedure we start 

9 The Ramanujan summation of the infinite series of even and odd integers: 

2- s + 4- a + ■■■= 2- s C(s), l- a + 3- 3 + ■ ■ ■ = (1 - 2~ s )C(s) 

leads to —1/6 and 1/12 for s = — 1, respectively [22] , 

10 The fact that P U {1} = P' U {2} implies that any successful regularization that can obtained with P can also 
apply to P', and viceversa. 



s 



with the two-point function of the field, given by the sum in modes 



1 00 1 

{0 R \<f>(t, X)$(t', sOlOfl) = [ e -in[(t-f-ie)-(x-x )]/R + e - in [(t-t' -ie)+(x-x )]/«] ? (39) 

2tt R 2vD' - 

n=l 



where, as usual, the infinitesimal ie has been appropriately inserted to make the two-point function 
well-defined. As is customary in quantizing a field in two-dimensions, we have discarded the zero 
mode n = in the above summation. At coincidence (t',x') —> (t,x) we have 



1 00 1 

(0 R \&(t,x)\0 R ) = -Y,-e- en/R 
Ztt * — ■ n 



(40) 



n=l 



which turns out to be the Abel summation of a divergent series [22]. From the sum 

00 

Y n- x e- en l R = - log(l - e~ t/R ) , (41) 



n=l 



one can easily derive the expression for the two-point function 

(0 R \$(t,x)$(j/,x')\0 R ) = -^-logfl - e -'(^-«)/«)(l _ e -i(Av-ie)/R, 

47T 



(42) 



where it is understood that the limit e — > + should be taken and An = u — u' := (t — x) — (t 1 — x') 
and Av = v — v' := (t + x) — (t' + x'). 

According to the above discussion, the energy density of the field, which classically is given by 



at J \ox J 
can be formally evaluated, at the quantum level, as 

oo 

(0r\T u \0 r ) = lhn (2vri2)- 1 J2 w ne~ en/R . 
In the limit e — > + , the quantity 

oo 

{2TiR)- l Y W nd~ en,R 



n=l 



n=l 



(43) 



(44) 



(45) 



is divergent, reflecting the well-known divergent vacuum energy. However, as is usual in quantum 
field theory one can invoke renormalization to properly find a finite quantity for the physical vacuum 
energy density. To renormalize (0 r \Tu\0r) one can perform an inverse radius expansion in 



(0 R \T$\0 R ) := {27iRy 1 Y J W n e 



en/R 



1 



1 



n=l 



2vre 2 24vri? 2 



+ 0(R~ 4 ) , 



(46) 



and renormalize by subtracting, from the exact expression (I45p . the terms in the expansion (jl6 
that diverge when e — > + 



(0 R \T tt \0 R ), 



lim 



1 



2ttR 2 {e^ 2R - e -^ 2R ) 2 2vre 2 



24vri? 2 



(47) 



In the CFT language, this result is somewhat related to the calculation of the central charge c = 1 
of the scalar field 1 131. 
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3.2 Scalar field with even/odd integer numbers modes 

It is now illustrative to consider the restriction of the modes of the original field <!>, parameterized 
by n, to the set of even (or odd) integer numbers. In this case we have 

oo 

(0 R \Tt en) \0R) = (2irRy 1 £ ^ e " m/i? = i^-^W + ° (jR " 4) ' (48) 

n=even 

and, after renormalization, we get 

(0 R \Tt en) \0R)ren = --^- 2 . (49) 

Similarly, for the quantum field defined on the odd integer numbers modes one gets 

(Q R \T^ d) \V R ) ren = -^. (50) 

We reproduce in this way the results obtained in the previous section with the zeta function 
regularization. 

3.3 Fermi field with (odd) prime modes 

Let us now extend the above analysis to the fields ip 1 introduced in section 2. Since, by definition, 
the (chiral) field ip (for simplicity we omit the index j) has only (odd) prime number modes in its 
Fourier expansion, the corresponding two-point function is then given by 

(0 R \i;(t - x)^(t' - x')\0 R ) = ^ E e-^-*'--)-^ ')]/« . (51) 

pGP 

At coincidence, one should deal with the following Abel sum 

(o R \^ ) \o R ) = ^Y, e ~ tp/R > ( 52 ) 

pGP' 

which, as expected, is divergent when e — > + . Moreover, in the calculation of the vacuum energy 
density one faces directly the Abel sum of all (odd) prime numbers 

(0nK nme (€) \0r) = £ P e^' R ■ (53) 

pGP' 

In sharp contrast with the field $, the restriction to the prime numbers turns the above sums into 
rather complicated functions of e/R. 

Our next and crucial step is to carry out a renormalization procedure to properly associate finite 
quantities to the above divergent sums. As in the previous analysis of the field this requires 
the introduction of suitable asymptotic expansions for large radius. The leading terms in the 
asymptotic expansions of (O^lT^^Oij) is essentially fixed by the corresponding ones in Minkowski 
space: (0 R \T$\0 R ) ~ l/27re 2 + ... . 

However, the definition of ip requires of a finite radius R and thus it is not possible to define ip 
in Minkowski space. Therefore, the asymptotic expansions needed to identify the subtracting terms 
cannot be based on the results in Minkowski space and one has to construct them according to new 
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criteria that we investigate in the following. Since our considerations should necessarily involve 
some conventional but sophisticated mathematics we will study the above Abel sums on the more 
conventional set P, instead of P'. Inasmuch as the fundamental aspects of the renormalization 
process only involves the behavior of the Abel sum as p — > oo, the shift from P' to P is irrelevant 
(one should add the quantity 1/A-kR 2 to the result obtained with P). 

Let us now focus on the renormalization of the expectation value of the energy density. To per- 
form renormalization we should expand (0_R|T^™ me ^|0r) for large radius. The desired expansion 
can be computed either directly or by differentiating the expansion of ^ P eP p e ~ ap ( a := e /-^) with 
respect to the parameter a twice and taking into account the additional factor — 1/ 4irR 2 . The result 
turns out to be the same. After an involved calculation it is possible to show that the asymptotic 
expansion of ^ P eP p e ~ aP w hen a — > satisfies (without invoking the Riemann hypothesis) 

log(-loga) + Bi + o(l) , (54) 

pGP P 

where B\ is Mertens' constant. This result is heuristically suggested by Mertens' theorem, but a 
rigorous proof involves non-trivial computations. For more details of this and other calculations 
we refer the reader to the Appendix at the end of the paper. 

Note that the growth of the first (and only) divergent term on the r.h.s. of Eq. (154p is very 
slow, due to the natural decrease of the density of prime numbers as p — > oo. To go further in the 
analysis we are forced to know the subleading terms in the expansion (154ft , For instance, terms 
of order 0(1/ log a) in the expansion of X^peP \ e ~ ap are converted, after derivation, into divergent 
terms. As explained in the Appendix one can obtain the following expansion if the RH is assumed 
to be true 



pTf p ti k ("log«) fc 

where T^ k \z) stands for the fe-order derivative of the T function. For example, the first five terms 
in this expansion are 

£i e -»p „ io g( -io ga)+ i 3l - ri i^- ' r2 ;^\ 1 

f^P (-log a) 12 (-log a) 2 

4C(3) + 7 7r 2 + 2 7 3 1 

6^ (-loga)3 + "' [bb) 

The expansion for (Oj?|T^ r " ne ^|0^) can be expressed then as 



This series has an infinite number of divergent terms (as e — > + ) and their sum does not converge. 
Therefore, one cannot define at this point a proper subtracting term to cancel out the divergences 
and, hence, one could be tempted to conclude that renormalization is not possible. 

However, a more careful inspection shows that this situation could be conceivably improved. 
To this end let us go back to the origin of (J55]). As explained in the Appendix, by accepting the 
truth of the Riemann hypothesis it is possible to derive the following asymptotic expansion 
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l e ~ ap = X - log(l + (- log a) 2 ) + B 1 - Im f°° e itlo ^ a (v{-it) - dt + o{a l l 2 ~ s ) , (58) 

that, by combining a further expansion of the first term with the third, leads to (|55p . Similarly, 

assuming the Riemann hypothesis, the asymptotic expansion for (0_R|T^" me ^|0j?) can be cast in 
the form 

W?r e (e) |0*> := ^ 2>~^ = ^(F(e/R) + o((e/Ry^)) , (59) 

pGP 

where all divergences in ([57|) are captured by the single function F(a) = o(a~ 2 ) 

F(a) = -Im— / e itlogc T(2 - it)dt . (60) 
a 2 Jo 

It is important to notice, however, that the terms of the form o(o _3//2_e ) in (|59p are potentially 
divergent. One would need to know the form of the non- leading contributions in (|59p to see whether 
or not the asymptotic expansion can be expressed as a finite sum of hierarchical divergent terms. 
In order to do this it would be necessary to go beyond the asymptotic analysis presented in the 
Appendix and study in detail the contributions to the asymptotic expansions coming from all 
the singularities in the integrands of (|62p and (|68p . The very complicated analytic structure of 
the relevant functions precludes us from making a definitive statement but, the presence of an 
infinite number of branch points originating in the non trivial zeros of the zeta function, strongly 
suggests that an infinite number of counterterms would be necessary to subtract all the infinities 
appearing in the expectation value of the stress-energy tensor. This would mean that the model is 
non-r enor malizable . 



4 Conclusions and final comments 

Motivated by the Goldbach and Polignac conjectures in Number Theory and the classical results on 
abelian bosonization in two-dimensions we study in this paper the "factorization" of a real scalar 
field (whose Fourier expansion is made out of even modes) as a product of two (so-called prime) 
fermion fields whose Fourier expansion exclusively contains prime modes. We have subsequently 
investigated the quantum properties of this theory by studying the vacuum energy density. As 
expected, the short-distance behavior of the prime field is highly non-trivial and quite different 
from a conventional field theory on the cylinder. Due to the restriction to the prime modes, the 
short-distance singularity is no longer purely of the form 1/e, but rather of form l/(eloge) + ... 
This feature makes the calculation of the renormalized expectation value of the energy density 
quite involved. The computation has turned out to be very elusive since an infinite number of 
divergences emerges naturally, showing the non-renormalizability of the underlying QFT. 

Let us stress that the Riemann hypothesis plays an important role in the calculation of the 
vacuum energy density. This is, in our opinion, a remarkable outcome of this work. The possibility 
of relating consistency properties of a particular field theory (such as the one considered here) to 
the Riemann hypothesis would definitely be very intriguing. 

The non-renormalizability of the theory -in the standard sense- would suggest that there is some 
scale at which its properties change dramatically. For instance, the two-particle sector of the Fock 
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space reflects the GC at low energies. Accordingly, the set of all two-particle states span all energies 
of the form n even /R. The non-renormalizability of the theory could signal a breakdown of the GC 
above some number density scale Nj R, acting as a sort of ultraviolet cut-off. By breakdown of the 
GC we should understand in this context that infinitely many even numbers cannot be written as 
the sum of two primes. 

On the other hand, since the value of the radius R can be freely chosen, the ultraviolet cut-off 
-and therefore N— should be in principle also as large as desired. This view actually suggests that 
the GC might be satisfied for arbitrarily large even numbers. One can envisage this possibility on 
deeper physical grounds by extending our study beyond QFT, e.g. to string-like models, where 
dualities involving the size of R (like T-duality, see for example [27)) could allow one to exploit this 
fascinating idea. 

Let us finally stress that no actual claim is made in this paper regarding the proof or disproof 
of the Goldbach conjecture. Our main result is to point out that a simple theory, a free field 
on a cylindrical spacetime with modes restricted to prime integers, is full of divergences, making 
impossible to properly define the vacuum energy. We believe that the truth of the Goldbach 
conjecture may be deeply related with this result, but in a way not unraveled yet. 



5 Appendix: Asymptotic expansions 

We give here some basic results about the asymptotic behavior as a — > + of the functions 

°° -ap 

/(»> = £ — 

oo 

g{a) = Y J PZ- aP 

where the sums are extended to the set of prime numbers P. As a way to highlight the role of the 
Riemann conjecture in the final form of the asymptotic expansions we will use complex variable 
methods. 



-ap 



5.1 Asymptotic expansion for f(a) 
The functional series 

oo 

V 

is convergent for a £ C with Re(a) > and uniformly convergent in every compact subset of this 
set, hence it defines a holomorphic function in the complex half plane {a S C : Re(a) > 0}. The 
best way to study the asymptotic behavior of / as a — > + is to use Mellin transform methods. In 
the present case the integral defining the Mellin transform converges if Re(s) > and is given by 

M[f; s] = V -M[e~ ap ; s] = T(s)P(s + 1) (61) 

per 

where P denotes the prime zeta function. This function is defined, in analogy with the standard 
zeta function, by the functional series 

1 



P(s) :=£- 



peP ' 
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for complex values of s with Re(s) > 1. The Mellin inversion formula provides us with the following 
integral representation for /. 



i pc+ioo 

f(a) = — a~ s T(s)P(s + l)ds (62) 

ziri Jc—ioo 

with c > 0. This representation is the main tool to obtain the sought for asymptotic formulas for 
/ and g. 

By using the expansion (see for example |23j ) 



P(s) = Y^^±\o E C{ks), (63) 



fc=i 



that can be obtained by using Mobius inversio we get 

1 i-c+ioo 00 ..(U\ 

f{a) = 2ri ■ a " T(s) 53 k bg C{k{s + 1)} ' 



that can be equivalently written as 



/(a) = — / a- s r(s)log(sas + l))ds-— / a"*r( S ) log( S )ds 

Z7T2 7 c—ioo ^ 7ri Jc—ioo 



' C—IOO 

+ ^-E / V/ a- s r( S )logC(Ks + l))ds. (64) 

Z7rz fc=2 K Jc-ioo 

for an appropriately chosen value of c. 

Let us start by considering the first integral in (j64|) . Without using the Riemann hypothesis 
this integral can be easily seen to be 0(1) by moving the integration contour to the line Re(s) = 0. 
This is possible because there are no zeros of the £ function on the line Re(s) = 1. Also, we have 
r(c + it) = 0(t c_1//2 e _7r * //2 ) and log|,s£(s + 1)| = 0(log|t|) and, hence, the contributions of the 
segments needed to close the integration contour and displace it to the line Re(s) = go to zero 
as L — > oo. The previous result can be improved by noting that the Riemann-Lebesgue lemma 
actually shows that the limit of the integral when a — > vanishes [or, equivalently, the integral is 
o(l)]. If the Riemann hypothesis is true the integral can be shown to be 0{a l / 2 ~ e ) for all e > by 
moving the contour to the line Re(s) = — 1/2 + e. Again this result can be slightly improved and 
the integral shown to be o(a 1 / 2 ~ £ ) as a consequence of the Riemann-Lebesgue lemma. 

Let us consider now the integrals appearing in 

— /-c+ioo 



1 n(h\ fc+ioo 

^E^/ a-°T{s)\o g ak{s + l))ds. 

k—2 Jc-ioo 



In each of the integrands, the singularities with real parts largei^l than -1 have several different 
origins: 



11 Here fx(k) denotes the Mobius function and £ is the standard zeta function. 

12 The only ones relevant for the present analysis. 
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• The ones with the largest real part correspond to the pole of T at s = 0. The residues are 
log CO) for k = 2,3, .... 

• The pole of £ gives rise to branch points at s = —1 + 1/k, with k = 2, 3, 

• The non-trivial zeros of the £ function, if the Riemann hypothesis is true, give rise to branch 
points with Re(s) = — 1 + 4-, and k = 2, 3, . . . 

As in the integral considered in the first step, the integration contour can be moved parallel to 
itself due to the exponential fall-off of the integrand and the mild growth of the absolute value of 
the C function. By moving the contour past the origin we get contributions, from the residues at 
the poles of the T(s) factor, that add up to the constant 

C:= f^m logC{kh 

k=2 

This, combined with the Euler constant 7, will give a final contribution equal to the Mertens 
constant B\ = 0.26149721 ••• . Once past the s = singularity, the integration contour can be 
displaced to the line c = — 1/2 + e and it is straightforward to show that the sum of the integrals 
on this contour is o(a 1 ^ 2 ~ <E ) without invoking the Riemann hypothesis. 



Finally let us consider the integral 

1 

~ 2vri 



c+ioo 

cTT(s) log s ds . (65) 



The singularities in the integrand are determined now by the poles of the gamma function at 
the non-positive integers, and the branching point of the logarithm at s = 0. The quick fall-off 
of the integrand provided by the gamma function allows us to freely move the integration contour 
in the negative direction of the real axis while wrapping around the singularities of the integrand. 
In order to compute the integral it is useful to somehow disentangle the singularity structure by 
choosing a branch cut for the logarithm such that arg(z) £ (— n/2, 3ir/2] (see figured]). As long 
as the choice does not change the value of the integrand in the region where the initial integration 
contour lies this is allowed. The advantage of proceeding in this way is that the poles of the gamma 
function will not lie on top of the branch cut. 

A straightforward computation gives now 

Re l"2^/ «" S r(s)log S d S j = 7+ -log(l + (-loga) 2 ) (66) 



Im / e in °z a (T(-it) - — )dt + 0(a 2 ) 
Jo * 



It is important to realize that despite the relative complication of (]66p it is actually easy to find 
its asymptotic behavior when a — > + . In particular, the last integral can be easily seen to be 
0(1/ (— logo)) by integration by parts and its full asymptotic expansion can be obtained directly 
by this procedure. This gives 



-Im I ^ { T H t) - — )dt ~ -Img ^—^^CTH*) - _) 

l"f. 

+ 



(_i)*r( fc >(i) 1 i^,(-i 



k=i v 6 1 k=i y & 
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Im(s) 




c + too 



Refsl 



c — too 



Figure 1: Integration contour for the integral (f65l) used in the study of f(a). 

In the limit a — > + the asymptotic behavior is controlled by the divergence of \ log(l + (— log a) 1 
that behaves as log(— log a) for small enough values of a. 



log(l + (- log a) 2 ) = log(- log a)-~Y, 



2^ k (-loga) 2fc ' 



By taking together the results obtained for each of the integrals in (|64p we conclude that, 
without invoking the Riemann hypothesis, we have 

/(o)=log(-loga) + Si + o(l), (67) 

whereas, if the Riemann hypothesis is true we have 

/(a) = -log(l + (-loga) 2 ) + J B 1 -Im / e H loga (r(-it) - — )dt + o{a^ £ ) 
2 Jo t 

for every e > 0. Here B\ := 7 + C is the Mertens constant mentioned above. As a side remark 
we would like to mention here that the asymptotic expansion (|67p is the same [11] as the one 
corresponding to the truncated series 

V -=log(-logo) + J3i + o(l) 
— ' p 

pGP;p<J- 



in the limit a — > 0. 
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5.2 Asymptotic expansion for g(a) 

The function g(a) can be studied by following the same steps as before. Now the functional series 

oo 

is convergent for a £ C with Re(a) > and uniformly convergent in every compact subset of this 
set, hence it defines a holomorphic function in the complex half plane {a € C : Re(a) > 0}. It is 
straightforward to find the following representation for g (a) by using Mellin transforms. 



1 



C+JOO 



9(a) = — I a- s T(s)P(s - l)ds (68) 



C — JOO 



with c > 2. This expression can be written as 

fc+ioo 



1 rc+ioo ^± (u\ 

g(a) = — a-T(.)X;^VlogC(fc(*-l))- 

Jc-ioo k=1 

The singularity structure of the integrand suggests to split the previous expression for g (a) as 
9(a) = — a- s T(s)log((s-2)as-l))ds-— a"T( S ) log(s - 2)ds 

ZlTl Jc—joo Z7TI Jc—ioo 

+ ^-E^l J / a- s r( S )logC(fe(s-l))ds. (69) 

zm k=2 K Jc-ioo 

for a certain value of c. 

Without invoking the Riemann hypothesis the first integral in (|69p can be easily shown to be 
0(a~ 2 ) by moving the integration contour; actually, by invoking the Riemann-Lebesgue lemma it 
can be seen to be o(a~ 2 ). If the Riemann hypothesis is true then this integral can be seen to be 
0(cr 3 / 2 - e ) (a bound that can, again, be slightly improved to o(a 3 / 2 e ) by using the Riemann- 
Lebesgue lemma). 

Let us consider now the sum of integrals 

f c+ioo 



1 ii (k) r c + l °° 

-E^V / a-T( S )logC(Ms-l))ds. 



2tri 

k=2 

At variance with the analysis presented for /(a), in this case the singularities of T(s) do not play 
any role to determine the asymptotic behavior of g(a). This expression can be seen to be o(a -3 / 2-£ ) 
regardless of the truth of the Riemann hypothesis. 
Finally let us study the behavior of 

1 /-c+ioo 

-— / a- s r( S )log( S -2)ds. (70) 

The singularities of the integrand and the the integration contours used to study the asymptotic 
behavior of this integral are shown in Fig. [5J A simple computation gives now 

(i rc+ioo \ i roo 
/ a- s r(s)log(s-2)ds = — 7 Im / e itloga T(2 - it)dt + 0(1) . 
27T« 7 c -ioo / « 2 Jo 
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The asymptotic behavior of this last integral can be easily obtained by integration by parts 



1 C OO "1 Jsi I " — ' 1 

~Jm / e itloga r(2-it)dt ~ -^ImV 
a Jo a k=0 



(-loga) fc + lL i { ' a? ^ (-loga)^ 1 ' 

fc=0 



Gathering the preceding results we conclude that without invoking the Riemann hypothesis we 
can only state that 

g (a) = o{aT 2 ) , 

whereas, if the Riemann hypothesis is true, we have that 

1 f'°° i 
g(a) = -Im-2 / e in ° sa T(2 - it)dt + o(a"2- £ ) 
a Jo 

for all e > 0. 

It is interesting to point out here that the leading asymptotics of g(a) can be obtained by 
differentiating twice the asymptotic behavior of the leading terms in f(a), i.e. 

1 f'°° ip~ l 

-log(l + (-loga) 2 ) + Bi-Im / e itloga (T(-it) - —)dt ~ 



log( - loga) + Bl -f izll^W / 

' fc (— logar 
fc=i 




-1 




Im(s) 




C + 100 



Re(s) 



Figure 2: Integration contour for the integral (|70p used in the study of g(a) 
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